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Abstract 

Large  Eddy  Simulations  (LES)  were  performed  to  study  the  flow  physics  and 
heat  transfer  in  a rotating  ribbed  duct.  This  numerical  study  is  simulating  the 
experiments  conducted  to  investigate  the  effects  of  buoyancy  and  Coriolis  forces 
on  heat  transfer  in  a turbine  blade  internal  coolant  passages  (Wagner  et  al., 
1992).  The  computations  were  performed  at  a Reynolds  number  (Re)  of  12,500 
based  on  average  velocity  in  the  duct  and  the  hydraulic  diameter  of  the  square 
duct.  The  rotation  number  ( Ro ) is  0.12  and  the  inlet  coolant-to-wall  density  ratio 
(Ap/p)  is  0.13.  The  rib  height-to-hydraulic  diameter  ratio  ( e/D ) is  0.1  and  the  rib 
pitch-to-height  ratio  (P/e)  is  10.  The  ribs  are  square  in  cross-section  and  are 
placed  transverse  to  the  flow  in  the  duct.  A direct  method  of  computing  a source 
term  due  to  unsteady  calculation  of  uniform  wall  temperature  case  in  periodic 
geometries  is  presented.  The  details  of  flow  field  and  the  temperature  field  are 
presented  and  analyzed. 

1.  Introduction 

Modern  gas  turbines  operate  at  very  high  turbine  inlet  temperatures  for 
better  second  law  efficiency  and  specific  thrust.  However,  such  increased  thermal 
loads  can  deteriorate  the  blade  life  in  a rotating  environment.  These  blades  are 
internally  cooled  by  using  the  serpentine  channels  with  turbulators  inside  the 
blade  to  enhance  the  heat  transfer  (Morris,  1981).  The  increment  in  heat  transfer 
due  to  rib  turbulators  as  compared  to  the  increased  pressure  drop  in  the  channel  is 
a crucial  design  parameter.  The  problem  is  complicated  further  due  to  the 
interplay  of  Coriolis  forces  and  buoyancy  forces.  Several  experimental 
investigations  to  study  the  effect  of  centrifugal  buoyancy,  rotation  number  and 
Reynolds  number  have  been  performed  (Wagner  et  al.,  1992).  However,  a 
numerical  study  can  provide  much  more  detailed  information  on  flow  physics  as 
well  as  heat  transfer  in  such  situations.  The  secondary  flow  in  non-circular  duct 
without  rotation  is  generated  due  to  the  anisotropy  in  turbulent  stresses.  In 
rotating  ducts,  the  Coriolis  forces  give  rise  to  secondary  flow  as  well.  Again, 
buoyancy  forces  can  generate  secondary  flow  field  to  enhance  the  crossplane 
mixing.  Low-level  turbulence  modeling  using  two-equation  models  is  incapable 
of  capturing  essential  physics  due  to  isotropic  nature  of  modeled  turbulent 
stresses  (Hermanson  et  al.,  2001).  Second  moment  closure  has  some  promise  in 
that  direction  (Jang  et  al.,  2001).  To  understand  the  unsteady  dynamics  of  various 
flow  structures  on  the  heat  transfer,  it  is  inevitable  to  use  LES  (Murata  et  al., 
2000  and  2001). 
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2.  Governing  Equations 

The  non-dimensional  governing  equations  for  conservation  of  mass,  momentum 
and  energy  for  an  incompressible  Newtonian  fluid  in  LES  methodology  are  as 
follows 
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temperature  and  Tt  , T2  are  yet  undefined  ref^rencf temperatures.  The  mean 
pressure  gradient  in  flow  direction  is  dP/dz.  Therefore,  p is  the  periodic 
component  of  the  pressure  field.  8y  is  the  Kronecker  delta  tensor.  £yk  is  the 
alternating  tensor.  The  distance  vector  can  be  written  as  rt  = RmSi3  + xt , where 
Rm  is  the  mean  radius  of  the  periodic  module  from  the  rotation  axis.  The 
important  parameters  for  such  flows  are  Reynolds  number  (Re  =UmDl/v), 
rotation  ^number  2(7?o  3 jQ . and  centrifugal  buoyancy  number 
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this  study.  Dynamic  Mixed  Model  (DMM)  is  used  to  model  these  SGS  stress 
tensor  and  scalar  flux  vector  (Moin  et  al.,  1991,  Vreman  et  al.,  1994).  The  box 
filters  are  used  in  the  Germano  identity  for  the  calculation  of  dynamic  coefficient 
and  for  the  calculation  of  Leonard  stresses.  The  dynamic  coefficient  is  test 
filtered  to  avoid  numerical  instabilities.  Treatment  of  the  non-dimensional 
temperature  in  the  periodic  direction  needs  special  attention.  Patankar  et  al. 
(1977)  described  a method  to  solve  the  uniform  heat  flux  (UHF)  and  uniform 
wall  temperature  (UHT)  problems  in  the  ducts  with  periodic  cross-sections  for 
steady  situations.  Wang  and  Vanka  (1989)  presented  an  iterative  procedure  to 
calculate  X.  However,  this  parameter  can  be  calculated  directly  for  explicit 
schemes.  Most  of  the  simulations  were  performed  using  non-dimensionalization 
with  respect  to  friction  velocity  and  uniform  heat  flux  case.  As  it  will  be 
explained  later,  that  renders  the  sink  terms  in  momentum  and  energy  equations, 
i.e.  dP/dz  and  X,  constant.  In  experiments,  usually  the  mass  flow  rate  and  wall 
temperatures  are  control  parameters,  therefore  reference  velocity  should  be  the 
average  velocity  and  X is  no  more  constant. 

Boundary  conditions  for  non-dimensional  temperature  in  the  periodic 
direction  is  written  as 
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where  superscript  indicates  the  z-location  and  subscript  b denotes  the  bulk 
non-dimensional  temperature.  Differentiating  the  periodic  boundary  condition 
in  the  wall-normal  direction  we  get 
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This  is  equivalent  to  enforcing  periodicity  on  the  Nusselt  number  in  a periodic 
geometry. 


Uniform  Heat  Flux  (UHF)  case: 
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Clearly,  setting  U to  Th  will  render  the  denominator  as  a constant.  Moreover,  the 
independence  of  non-dimensional  bulk  temperature  from  periodic  direction 
implies  that  1)  is  equal  to  Th.  Therefore,  the  scaling  at  the  inlet  plane  and 
periodicity  of  Nusselt  number  can  uniquely  determine  the  non- 
dimensionalization  and  the  sink  term  in  the  energy  equation. 


Uniform  Wall  Temperature  (UWT)  case: 
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Here  q is  the  wall  normal  direction  and  dS  is  the  differential  area  element  on  the 
wall. 

For  the  square  channel,  we  use  Nusselt  number  periodicity  to  define  the  flux  at 
the  inlet  in  terms  of  the  flux  at  the  exit  as 
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Using  these  relations  in  the  energy  balance,  we  get 
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To  enforce  the  validity  of  scaling  relation  up  to  the  wall,  we  choose  Tt  equal  to 
Tw.  Therefore,  the  non-dimensional  temperature  is  zero  at  the  wall  and  the 
scaling  ensures  the  periodicity  of  the  Nusselt  number  in  the  periodic  geometries. 

3.  Problem  Description  and  Computational  Method 

The  computations  were  performed  at  a Reynolds  number  (Re)  of  12,500 
based  on  average  velocity  in  the  duct  and  the  hydraulic  diameter  of  the  square 
duct.  The  rotation  number  (Ro)  is  0.12  and  the  inlet  coolant-to-wall  density  ratio 
(Ap/p)  is  0.13.  The  rib  height-to-hydraulic  diameter  ratio  ( e/D ) is  0.1  and  the  rib 
pitch-to-height  ratio  (P/e)  is  10.  The  ribs  are  square  in  cross-section  and  are 
placed  transverse  to  the  flow  in  the  duct.  This  numerical  study  is  simulating  the 
experiments  conducted  to  study  the  effects  of  buoyancy  and  Coriolis  forces  on 
heat  transfer  in  turbine  blade  internal  coolant  passages  (Wagner  et  al.,  1992). 


Figure  1 Schematic  of  the  computational  domain 

4.  Results  and  Discussion 

The  three-dimensional  spectrum  of  the  instantaneous  flow  field  is  shown 
in  figure  2a.  The  grid  resolution  is  sufficient  to  capture  the  energy  producing 
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events  as  well  as  the  portion  of  the  inertial  subrange.  A peak  in  energy  spectrum 
is  also  observed  around  the  wave  number  corresponding  to  a length  scale  l/D  = 
0.16.  Clearly,  this  can  be  attributed  to  the  energy  production  by  vortex  shedding 
behind  the  ribs  ( e/D  = 0.1).  To  maintain  the  flow  rate,  the  mean  pressure  gradient 
is  applied.  However,  this  pressure  gradient  is  subjected  to  temporal  variations 
corresponding  to  the  Strouhal  frequency  of  vortex  shedding  due  to  ribs.  This 
pulsation  causes  the  flow  rate  to  vary  with  the  same  frequency  (figure  2b).  The 
average  flow  rate  is  maintained  at  1 .0  as  desired. 


Figure  2 (a)  Three  dimensional  energy  spectrum  of  the  flow  field  (b)  Flow  rate  vs 
time. 

FLOW  FIELD  DESCRIPTION 

The  instantaneous  snapshots  of  the  streamwise  component  of  vorticity 
field  at  three  different  streamwise  stations  (z/D  = 0.25,  0.50  and  0.75)  depict  the 
complex  flow  field  in  the  duct  (figure  3).  At  z/D  = 0.25,  the  rib  is  placed  on  the 
leading  wall  and  sheds  vortices  into  the  mainflow  with  significant  streamwise 
vorticity.  The  boundary  layer  on  the  trailing  wall  is  highly  turbulent  and 
intensified.  The  local  increase  in  velocity  magnitude  near  the  trailing  wall  is  due 
to  decrease  in  cross-sectional  area  as  well  as  body  forces.  At  the  streamwise 
centerplane  z/D  = 0.50,  the  recirculation  region  behind  the  rib  on  the  leading  wall 
contains  intense  vortices  along  the  flow.  The  trailing  wall  vortices  are  gathered 
towards  the  center  of  the  duct  due  to  secondary  flow  in  the  crossplane  of  the 
duct.  At  z/D  = 0.75,  the  rib  is  placed  on  the  trailing  wall  and  interacts  with  the 
oncoming  turbulent  boundary  layer  vortices.  In  this  snapshot,  the  vortices  are 
pushed  towards  the  center  of  duct  by  the  rib  as  well  as  the  secondary  flow.  The 
leading  wall  boundary  layer  shows  a lot  of  activity  too.  The  vortices  at  the  front 
and  the  back  wall  do  not  penetrate  into  the  core  flow  to  the  similar  extent  as  the 
leading  and  trailing  wall  vortices. 
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z/D  = 0.25  z/D  = 0.50  z/D  = 0.75 

Figure  3 Streamwise  component  of  instantaneous  vorticity  CD*  at  three  XY  planes 


Time-averaged  velocity  field  shows  the  skewed  profile  (Figure  4).  The 
boundary  layer  on  the  trailing  wall  (unstable)  is  much  steeper  than  on  the  leading 
wall  (stable).  The  details  near  the  ribs  show  the  difference  in  the  size  of 
recirculation  regions  in  front  and  behind  the  ribs.  The  flow  attaches  over  the  top 
of  the  rib  B (on  the  trailing  wall)  while  it  remains  detached  on  the  top  face  of  rib 
A (on  the  leading  edge).  The  impingement  of  the  oncoming  flow  on  the  front  of 
the  ribs  results  in  high  heat  transfer  rates.  The  recirculation  behind  the  ribs 
accompanied  with  enhanced  crossplane  mixing  will  result  in  the  increased  heat 
transfer  rates. 


Rib  A 


Time-averaged  Velocity  Vectors  at  Y/D  = 0.5  Rib  B 


Figure  4 Time-averaged  velocity  vectors  and  details  of  flow  field  near  the 
ribs  at  the  Y/D  = 0.5 
HEAT  TRANSFER  RESULTS 
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Instantaneous  non-dimensional  temperature  field  is  shown  at  the  first  cell 
node  over  the  corresponding  walls.  There  is  more  coolant  accumulation  on  the 
trailing  wall  as  compared  to  the  leading  wall  (note  the  difference  in  the  range). 
The  heat  transfer  is  enhanced  on  the  trailing  wall  by  a factor  of  two 
approximately.  The  temperature  field  show  streaks  correlated  with  the 
streamwise  component  of  vorticity  on  these  walls.  The  coolant  fluid  in  front  of 
the  ribs  increases  heat  transfer  in  the  stagnation  (front  recirculation)  region.  In 
the  leeward  recirculation  regions,  the  non-dimensional  temperature  is  close  to 
wall  temperature.  The  temperature  distribution  on  the  front  and  back  wall  is 
similar.  However,  there  is  more  coolant  near  the  back  wall  as  compared  to  front 
wall  (and  it  is  observed  through  out  the  computational  duration).  This  might  be 
caused  by  a low  frequency  mode  in  the  coreflow.  Again,  the  temperature  field 
correlates  with  the  streamwise  streaks  of  coolant  fluid  on  these  walls. 
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a)  Trailing  wall 
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b)  Leading  wall 
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c)  Back  wall 


d)  Front  wall 


Figure  5 Instantaneous  temperature  field  on  the  walls  of  the  duct  a)  Trailing  wall, 
b)  Leading  wall,  c)  Back  wall  and  d)  Front  wall. 

5.  Conclusion 

Large  eddy  simulations  were  performed  for  a rotating  square  duct  with 
normal  rib  turbulators  to  enhance  the  heat  transfer.  The  Coriolis  force  as  well  as 
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centrifugal  buoyancy  parameter  has  been  included  in  this  study.  A direct 
approach  is  presented  for  the  unsteady  calculation  of  non-dimensional 
temperature  field  in  periodic  domains.  The  complex  flow  field  shows  dominant 
secondary  flow  vortices  that  enhance  the  mixing  of  the  thermal  boundary  layers 
on  the  duct  walls  with  the  coolant  fluid  in  the  core.  The  temperature  field  is 
highly  unsteady  and  may  contain  low  frequency  mode  that  allows  the  coolant  to 
adhere  to  either  front  or  the  back  wall  of  the  duct.  In  the  future,  the  assessment  of 
heat  transfer  enhancement  with  respect  to  increased  pressure  loss  will  be 
presented. 
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